Boson Interferometry after SPS and before RHIC 
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We re-examine the connection between interferometry and the Wigner representation for source 
freeze-out, and continuous emission. At the operator level, two equivalent representations of the 
two-particle spectrum are found, which contradict the standard expression of kinetic theory. The 
discrepancy is resolved using two toy models. Further, we revisit interferometry in scale-invariant 
one-dimensional hydrodynamics, and argue that recent experimental results are evidence for a short 
kaon emission time. Using two exactly calculable models of two- and three-dimensional flow, it is 
shown that the saddle point approximation, which is reasonable for one-dimensional flow, is no 
longer adequate. In these models the scaling law is altered, and we argue that such qualitative 
trends, together with other observables, are vital if one is to draw conclusions about the unknown 
source parameters. 



1. Introduction 



The hope of discovering a quark-gluon plasma (QGP) in heavy ion collisions is to some extent connected to 
possibility of measuring the geometric size of the region of secondary particle production. A quantitative estimate of 
this size is necessary to obtain the energy density, an important quantity in the discussion of the deconfinement phase 
transition. An important tool in accomplishing such a size measurement is interferometry. 
, Recent experiments at the CERN SPS indicate that interferometry is sensitive to hydrodynamic motion in the 
' emitting source At RHIC energies one may expect the signature of this motion to be even more pronounced. 

T-H , In this paper, we therefore wish to reexamine various techniques to calculate the two-boson spectra for expanding 
Tij" • sources. 

I Intensity interferometry was proposed by Hanbury-Brown and Twiss to measure stellar sizes. The nature of light 
I emission in stars is well understood; the thermal mechanism guarantees that photons are emitted independently from 
^ i' different parts of the photosphere, while the small angular size of the star makes it impossible to obtain an optical 
' image. These two factors create the conditions necessary for interference between the two-particle amplitudes of 
photons emitted by different parts of the stellar surface. Under certain assumptions which will be discussed later, the 
(~| , two-particle detection probability is 



W{ki,k2) = / d'^xid^X2 pixi) p{x2) 



1 + cos(fci ~ k2)-{xi - X2) 



(1.1) 



X 

I The nature of particle sources in nucleus-nucleus collisions [p| is less well understood. In particular, the correspon- 



dence between measured quantities and parameters of the emitting system is less clear. One possible reason Eq. (1.1) 
may be inapplicable is the presence of correlations on the same scale as the size of the emitting system. An extreme 
case - when interferometry is completely unrelated to the total source size - is the scattering of coherent light in 
fluids; the rate of two-photon coincidence depends on the three-particle distribution function in the fluid Q]. For 
A-A collisions, even in the simplest scenario of hydrodynamic evolution, the relation between the inclusive one- and 
two-particle spectra and the parameters of the emitting system does not follow the classical scheme of interferom- 
etry. The width of correlator is not directly connected to the geometric size of the source, an effect first discussed 
in Ref. For an expanding spherical shell, the apparent (or "visible") source size is smaller, the larger the total 
pion pair momentum P±. In nuclear collision models, the dependence of the apparent size on P± may reflect a more 
complicated interplay between the collision dynamics and the true source size |^ . 

A relativistically covariant theory for interferometry in the hydrodynamic model of nuclear collisions was developed 
in Ref. JR] , with the main focus on freeze-out as a realistic mechanism for final state hadron production. In a subsequent 
study [pT, it was shown that if the dependence of the correlation function on the difference of the longitudinal momenta 
is rescaled by y'm±, the width of the correlator becomes independent of m^. This scaling behavior has also been 
discovered in parallel studies |^ based on the intra- nuclear cascade approach jl^, where classical currents are the 
sources of final-state pions. The m± scaling dependence has been observed at SPS energies by the NA35 and NA44 
collaborations [lli 
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A covariant theory of interference for continuous emission has also been applied to dilepton and photon [g2 13 
emission from a quark-gluon plasma. Analysis of polarization effects in photon interference |12P indicates that there 
is a dilemma in the proper choice of formalism: the languages of locally defined states, and Wigner phase-space 
distributions lead to different answers. 

Events generators such as RQMD [Q suggest a classical description of the space-time evolution of high energy 
heavy ion collisions. Other models, based on the solution of semi-classical kinetic equations, often provide descriptions 
in terms of one-particle distributions. Interferomctry is a manifestly quantum phenomenon, and we shall therefore 
begin in Section 2 by examining its connection to the widely used Wigner representation, for the cases of source 
freeze-out (initial data problem) and photons (continuous emission). At the operator level, the Wigner representation 
is simply a formal re-expression of the one- and two-particle spectra, and we obtain two equivalent expressions for 
the two-particle spectrum. These are in contradiction with the standard expressions of kinetic theory. In Section 3, 
we attempt to resolve the issue by way of two toy models, viz.^ the emission of particles from one and two cavities. 
In an appendix, we reanalyze the derivation of the two-particle spectrum in classical source models. 

In Section 4, we shall revisit interferometry in scale-invariant one-dimensional hydrodynamics, deriving the mj_- 
dependence of the longitudinal effective size in the case of the "extended freeze-out." We argue that recent experi- 
mental results confirming this behavior are evidence for short kaon emission times, and in disagreement with the 
prediction of kaon distillation . 

Unlike the inverse scattering problem, where an analysis of the phase shifts allows one to obtain the shape of the 
potential, the inverse problem of interferometry has no exact mathematical formulation. The interpretation of the 
two-particle spectra requires a model that is determined a priori, up to the numerical values of its parameters. We 
therefore present in Section 5 some results for two- and three-dimensional flow in exactly calculable examples: the 
transverse explosion of a long filament, and the spherical explosion of a point-like source. In both cases the saddle 
point approximation, which is reasonable for one-dimensional flow, is no longer adequate. The scaling law is altered, 
and we argue that such qualitative trends, together with other observables, are important in establishing the type 
of model one is dealing with. Only once the model is given, can one draw conclusions about the unknown source 
parameters. We conclude in Section 6. 



2. Initial data and continuous emission for interferometry 

The collision dynamics in heavy ion physics is rather complicated. However, some stages of the collision may be 
described in terms of one-particle distributions, or even by a few macroscopic parameters. A strong, yet attractive, 
simplification occurs if we assume that the hot expanding matter freezes out at some critical temperature Tc- Generally, 
such models rely on (semi-) classical assumptions, and produce (semi-) classical final distributions. On the other 
hand, interferometry is a manifestation of the quantum nature of the constituents of the system, and it is therefore 
necessary to examine carefully whether these distributions can serve as an adequate input. 

A natural connection between quantum mechanics and (semi-) classical descriptions is provided by the density 
matrix p. In a formal language, the problem we wish to solve in interferometry is to determine judiciously chosen 
parameters of p by measuring the inclusive cross-sections dNi/dk, dN2/ dkidk2, .... 

2.1 Source freeze-out 

Let |in) be one of the possible initial states of the system emitting a pion field '^{x). At this stage, we do not wish 
to specify the exact nature of the initial state, but study instead the observables relevant for interferometry at the 
operator level. 

The field 'f>{x) is to be detected later by some device, or analyzer. Let it be tuned to the measurement of the 
momentum of the free pion. Thus, the eigenfunctions of the analyzer are the free pion wave functions 

f^{x) = (2^)-3/2(2fcn)-i/2e-"=- . (2.1) 

These functions are defined for > tc, where tc is the time when the system "prepared" the final pion. The 
corresponding annihilation operator Aj: for momentum k in the final state is given by 

4- = I d^xJl{x)idlm , (2.2) 

where (ad^b) = o-id^b) — {d^ a)h. The operator describes the effect of a detector (analyzer) far from the point 
of emission, so, by definition, the pion is detected on mass-shell, = (J? -\- vn?)^/'^. The inclusive amplitude to find 
one pion with momentum k in the final state is 



2 



{X\Aj:S\m) , 



(2.3) 



where S is the evolution operator after freeze-out, and the states \X) form a complete set of all possible secondaries. 
Summing the squared modulus of this amplitude over all (undetected) states and averaging over the initial 
ensemble, we find the one-particle inclusive spectrum 



dNi 
dk 



(2.4) 



where the density operator will describe the emitting system. For now, we shall study only the operator part of 
Eq. (2^), N-^ — At^^, which gives the number of pions detected by an analyzer tuned to momentum k. In the same 
way one may obtain the inclusive two-pion spectrum 

dN2 



= Trpin^t4 4 AtAt S 
dkidk2 ' ' 

and the "device" operator which counts the number of pions pairs is 

it 



(2.5) 



(2.6) 



The pion field <p(x) reaches the detector after free propagation. This simplest kind of evolution is described by the 
retarded Green function 

= J dT.f,{y) G,.ot(x - y) dy 0{v), 

^\x) = J dS^(y) ^t(y) G^^Ay - ^) , (2.7) 

where the pion field ip{y) is given on the 3-dimensional freeze-out hypersurface (a Cauchy hypersurface) . The space 
of states in which the density matrix acts is also defined on this surface. Substituting Eqs. (2.7) in (2.2), we find the 
final form of the pion Fock operator, expressed in terms of the initial fields: 



S^Aj^S = j d'xflix)id, J dJ:^iy) Grotix ^ y) dy viv) , 
S^AlS^ f d^x [ dJ:^{y) ^t(y) Gadv(2/- x)*a! /fc(x) 



These equations may be simplified using the explicit form of the free pion propagators: 



Grct {x - y) 



d^p 



-ip-{x-y) 



(2.8) 



(2.9) 



S^A,S = e{x° - y°) / dl]^(y) my)il)y 0{y) , 



(27r)4 (pO ± ie)2 - p2 _ ^2 • 
Substituting Eqs. (^.9|) into (^.8D and integrating, first over d^x which sets p equal to k, and then over we obtain 



(2.10) 



As before, the pion four-momentum k is on mass-shell. This is not surprising, as the propagation is free and our 
analyzers perform an on -shell Fourier expansion of the initial data. The 0-function appears because the propagation 
is retarded. Using ( 2.10 ), we may rewrite the number operators for single pions and pairs of pions as 



S^A\S = 0{x" - y") / d^^iy) ^\y)i^. My) 



N. = eix"-y^) d^Md^M 



fk{yi)idy, <^Hyi) 



<f>{V2)idy^ fl{y2) 



(2.11) 



and 
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(2.12) 



respectively. 

We shall now attempt to incorporate the technique used in kinetic studies by passing to the Wigner representation 
at the operator level, i.e., without specifying any physical information about emitting system. Firstly, for Nj^, we 
choose new variables 



2/1 = ^+ 2 ' y2 = R - ^ 



(2.13) 



This transformation is motivated by a simple idea: we expect that R will label the sources, while the Fourier transform 
over z will yield the local spectrum. We assume that the 3-vectors R and z lie in the Cauchy hypersurface, which is 
obviously not planar, in general. The local time-like direction coincides with the normal. As the equation for ip is 

^0 

second order, the initial data includes the derivatives dy., which lie along the local normal. To adjust the variables, 
we write 



dyi ^ -^dB. + dz , 



^0 1^0 ^ 
9y2 = 2^^^ ^ 





z ' 



(2.14) 



allowing us to introduce a compact form of the Wigner representation of the operator (2.11) 

^0 



l3., „~i(ki+k2)-z/2 ~ik°-{F!°+z°/2) 



X(^(i?"-i-,i?-^)* 



1 ^0 ^0 

2 dn + dz 
2 dR- 



^ni?- + -,i?+-) 

^^k°■iR'>+z'>/2) 



(2.15) 



which may be shown to depend only on the combination (fc5' + fc2)/2. If the Cauchy surface is planar, all the derivatives 
above will be along the common normal. Then, the one particle spectrum may be represented in a compact symbolic 
form 



Nk 



1 
2fcO 



d^RAf{R,k) 



(2.16) 



For the two-particle operator we have four arguments in the integrand of (2.12), which allows us to choose the 
variables for the Wigner transformation in two ways: 



^ ^ zi ^ Zl 

yi ^ Ri + —, t/2 = -Ki - —, 

2/3 = «2 + y, 2/4 = -K2 - y 



(2.17) 



or, equally. 



yi = ^1 + y, 2/4 = ^1 

2/3 = ^2 + y, m = R2 



2 ' 

f2 

2 



(2.18) 



Correspondingly, the operators for the two-particle spectrum may be rewritten in the two identical forms 



1 



--^l^i -^2)) = -^^^\d^Rid^R2M{RiM)M{R2M) 



(2.19) 



and 



d^R.d^R, M{R^, h±^)M{R^, ^1±^) e-^Ck.-k.)iR^-B.2) 



(2.20) 
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respectively. 

2.2 Continuous emission 

We shall now consider continuous emission from an evolving system, using the case of photons as an example . 
Unlike for pion emission, where the (transition) currents which might emit the pions are not conserved, this problem 
has a clearer footing. Also, technically, the Wigner representation for continuous emission is simpler, as there are no 
derivatives in the equations that will replace ( ^.11 ) and (2.12). As in the previous section, we define the inclusive 
probabilities for one and two photons 



dN- 



dk 



and 



dN-, 



Trpin 5''^c''^(A;i,Ai)c'f(fc2,A2)c(A:2,A2)c(fci,Ai)5' , 



dki dko 1 N 

respectively. The operators c{k) appear in the decomposition of the electromagnetic field operator. 



(2.21) 



(2.22) 



(2^)3/2 (2c^,-) 1/2 



c(/c,A)e-*'^-^ +ct(fc,A)e*''-=^ 



(2.23) 



where A runs over the two physical polarizations of the photon. Since the quantum equation of motion for the field 
A^(x) has the classical form, 



A^(x) = d'yDZ{x-y)My) , 



it is easy to show that the operator of the number of emitted photons is 

d'^xd'^y 

X 



2fcO(27r)3 



-ik-{x-y) '-t 



jl{x)iv{y) 



Now we may safely perform a full 4-dimensional Wigner transformation 



(2.24) 



(2.25) 



Introducing 



we may rewrite 



:-t 



2fcO J {2tt 



3l{R + z/2)UR-z/2) 



Y,el'\k)ei'\k) f d'Rn'^-'{R,k) 



For the operator that counts the number of photon pairs, 



The two-photon detection function can be rewritten two ways, by changing variables either as 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



(2.30) 
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^1 r, 



2 2 ^2 



(2.31) 



These transformations lead to two equivalent representations of the two-photon spectrum: 



(2.32) 



and 



N^,^iN^,^~S{k,-k2)) = 5] e(^^)(fci)6(^^)(fci)4'^Hfc2)4'^nfc2) 

d^Rld^R2 m^{Rl, h±h) nP-(i?2, .^iLl^) COs(i?i-i?2)-(fci-fc2) 



respectively. 
2.3 Discussion 



(2.33) 



We have obtained three different expressions for the two-particle spectra, for both the initial data (freeze-out) 
problem and for continuous emission, viz, Eqs. ( p. 12 ), (2.19) and (2.20), and (2.29), (2.32) and (2.33), respectively. 
An apparent paradox appears if one compares these expressions with those commonly used. For the freeze-out 
mechanism of the pion production these expressions have the form 



4fc?fc°4^- / d^Rid'^R2N{Ri,ki)N{R2,k2) 



dkidk2 



(2.34) 



d^R^d^R2 AA(i?i, h±h) AA(4, h±hl) ^-^ik^-k.)■iR^-R.) 



while for continuous photon emission 



4fc?fc0-^= E ^^,''\k)ef^\k,)e^^^\k2)ei'^\k2) f d^R,d^R2 
dkidk2 , , J 



(2.35) 



Ai ,A2 

n^''(i?i,A;i)n'''^(i?2,fe) 



(For references and a detailed derivation of such formulae, see Ref. |l|l; there, the off- mass-shellness of the distribu- 
tions in the cross-term is emphasized as a physical result.) These expressions appear to be in line with kinetic theory: 
they seem to allow one to use the (semi-) classical phase-space distributions as the input f or su bsequ ent st udies of the 
interferometry problem. Below, we shall resolve the manifest contradiction between Eqs. (2.34) and (2.35) on the one 
hand, and the expressions given in the previous sections on the other. For now, we emphasize that the phase-space 
densities with argument (fci -I- fc2)/2 cannot be considered as distributions of the single particles in definite quantum 
states. For the photons, there appears an additional discrepancy between the transversality of the polarization tensor, 
and the physical polarization of the detected photons p2[ . 



3. Interferometry 

Up to this point, we have only discussed how the particles are detected, and have yet to describe the states in which 
particles were created. We have seen that the formal Wigner transformation puts the one-particle operators into 
a form resembling classical phase-space distributions, but that the representation for the two-particle "detection" 
operators is not uniquely defined. It allows for two equ ivale nt forms, each of which, by coincidence, resembles one of 
the terms of the familiar interferometric formula, Eq. (2.34), or Eq. (2.35). Each Wigner operator representation is 
supposed to yield the two term interferometric formula as a final result. 

To proceed with a discussion of interferometry, we require physical input. This input is already present in equations 
like (^^) and (^^), where the operators are averaged with the density matrix of the initial state. It is not evident 
that any arbitrary input will make interferometry which yields physically useful information possible. 
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Interference of quantum mechanical amplitudes appears when for identical initial states and identical final states 
(of one or more particles) there exists more than one intermediate transition amplitude. Interference in one-particle 
propagation exists, and is used in optical devices to obtain an optical image of the source. Interference in the 
propagation of the two-particle amplitude does not necessarily exist, because there may not be two alternative paths 
for the evolution of the two-particle state. An example is two-photon emission from a single atomic transition, and 
for each different case, we must examine the initial state of the system, i.e., the density matrix 



3.1 Toy model 1. Emission from a one-dimensional cavity 

Consider the case of an initial state which is a one-dimensional cavity filled with hot bosonic radiation. For the 
sake of simplicity, we shall impose periodic boundary conditions on the walls of the cavity. The cavity eigenfunctions 
are 



d,p{y) = (2a)-i/2(2/)-i/2e-'P !' , tor - a < y < a , 
= , otherwise 

where 2a is the distance between the walls, and the allowed values of the momentum p are 

p„ = — , n = 0,±1,±2, ... . 
a 



(3.1) 



(3.2) 



The boson annihilation and creation operators in the cavity are defined by their field decomposition: 

V{y) =^ap(j)p{y) , (3.3) 



so that, i f the walls are removed at time — 0, the average number of particles with the momentum k is given by 
(see Eq. ( pl|) ) 



dyidy2 ^n{p) 



My 2) dy^ /fe(2/2) 



(3.4) 



where n{p) — (ajop) is the boson occupation number. The integrals may be evaluated to yield 

{Nk) = > n{p) — — . 



Using Eq. (2.12), we may similarly write for the two-particle spectrum 

,,7 ..7 c x\ i^^k,+u)p,){u!k,+uJpJ{u;k^+u}p^){ojk2+uJpJ 



^ ('^pi '^P4 ^^^2 ) 



dye 



Pl ,--,P4 

i (ki-pi)y 



Air^Aa'^AuJki Wfcs V^Pi ^p3 '^P4 



dye 



i iki-p2)y 



dye 



-i (fe2-P3)y 



dye 



i (fe2-P4)y 



In case of the Gibbs ensemble, the average of the four Fock operators is easily calculated 



{al^al ttp^ap^) = n{pi)n{p3) 



Spip2^P3P4 + ^pipi^p 



SO that 



{NkANk2 -4ife2 = ^fei ^fc2 + 2^ n{pi)n{p3) 



Pl,P3 



sin(pi — ki)a sin(pi — k2)a sin(p3 — ki)a sin(p3 — k2)a 
{Pi-h) {Pi-h) {P3-ki) (^3-^2) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 
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Let us first discuss the single particle spectrum, (p.5|). From Eq. (3.2), it is clear that (Nk) has an infinite sequence 
of alternating zeros and peaks, in analogy to an optical Fabry-Perot interferometer, where the oscillations arise from 
the first order (one-particle) multi-ray interference |l^ . Here, we rely on the interference between all possible histories 
of a single particle, and, in principle, the distance between two neighboring zeros of (Nk) will determine the size of 
the cavity. 



The two-particle spectrum, Eq. (3.8), has a rather complex structure. It arises because any given two-particle 
final state can be traced back to an infinite number of initial states; when the cavity is opened, every eigenstate 
of the finite-sized cavity has, in addition to its "fundamental harmonic," an infinite sequence secondary peaks, or 
"satellites." A given plane wave particle may thus have originated from any single-particle level in the cavity. The 
resulting unwanted structure is an artifact of our choice of boundary conditions, and would dissapear if the zeros in 
(Nk) were not equidistant. 

As long as the detailed behavior inherent in this model is not expected to appear in a real physical applicati on o f 
interferometry, it is reasonable to modify it by eliminating the rapid oscillations in the two-particle spectrum Eq. ( [3.8[ ) . 
As artificial as the mod el is , any procedure will do. Ad hoc, one could remove the secondary spectral lines by hand, 
restricting the sum in ( [3.^ ) to \pi — kj\ < n/a. This scheme is essentially equivalent to taking a sufficiently large 
cavity (for a given momentum) - two of the four integrals in Eq.( [3.7| ) become delta-functions, and the other two can 
be written as 

sin^(fci ~ k2)a 

(fci - 

This factor is strongly peaked at ki = k2, with a corresponding suppression of the satellites; physically, for reasonable 
a, the detected particles can be traced back to their parent levels in the closed cavity, and no interference is possible. 
For fci = k2, we recover the two particle normalization. 



3.2 Toy model 2. Emission from two one- dimensional cavities 

Next, we consider two cavities, defined by walls at L ± a and —L ± a. There are now two sets of the eigenstates, 
defined separately for each cavity: 

0p,jv(2/) = (2a)-i/2(2fcg)-i/2e-»P-y ^ ioY Lm - a < y < Lm + a 

— , otherwise. (3.10) 



The index N runs over the two cavities, and Ln = ±L. The spectrum of states in each cavity is defined by Eq.( [3.5[ ), 
and the field decomposition is 

= ^o,p^N (l)p,N(y) , (3.11) 

p,N 

where the Fock operators have acquired an additional index enumerating the cavities. States belonging to different 
cavities are independent, as implied by the commutation relations: 

l&p.N, Op', A"] = ^pp'^NN' ■ (3-12) 

With the agreement that there is a one-to-one correspondence between the states in any one cavity and states in 
which we detect the particles, the one-particle distribution reads 

(Nk) = ^n(fc,Ar) - ^ (a[.,^a,,Ar) . (3.13) 

N N 



We now consider the two-particle spectrum, and our intention to resolve the discrepancy between Eqs. (2.12), (2.1£) 
and (^!20|) . In our model, there is interference between four amplitudes: two amplitudes when each particle is emitted 
from a different cavity, and two more when both particles originate from the same cavity. The general expression for 
the two-particle spectrum reads 

/jCr fjCr X \\ \^ /'t -t - - (^fc, 4- CJp J (^fc, -I- UJp^ ) (Cjfca +0-^3 ) (cJfc, + LJp, ) 

(iV..(A^..-4..J> = ^l^J%,,N.a^^.N.apM.,N.) 4.24^24^^^^^^^/,:,^== 

Ni[...\N4 

Li-\-a nL^^a nL^-^-a pL4-\-a 

dy e-'^'^^-P^^y dye'^^^-P^^y I dy e'^^'^^-P^^y dy e'^^^'^^^y , (3.14) 
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where Lj_= L{Nj) = ±L, and momentum pj originates from the cavity Nj. The statistical average may be found 
using ( 3.12 ): 



i"■l^,N^^U,N3^Pi,N4ap2,N2) = n{pi, Ni)n{p3, N3) 



Spip2^NlN2^P3P4^N3N4 + Spip4^NlN4,Sp3P2^N3N2 



(3.15) 



Eq. (3.14) becomes 

{K{Nk2~5k,k2)) = {K) {K) + 



5: n{p,,N,)n{p,^ m) ^""'^ +^P.)i^k. +^PJK. +-Ps) 



PI ,P3 



47r24a24u;fejWfc2Wpja;j 



P3 



^^ sin(pi - ki)a sin(pi - k2)a sin(p3 - ki)a sin(p3 ~ k2)a g-.,(Li-L3)(fei-fc3) 
{Pi - ki) {pi - fe) (P3 - fci) {P3 - fe) 



(3.16) 



The two terms in the sum with iVi ^ N3 {i.e., Li ^ L3) lead to the usual interference term, with periods of oscillation 
Afc ~ 1/(2L). 

The resolution of the detectors should not allow us to determine which cavity the particle originated from, i.e. we 
do not have 



\ki - fel i > 1 



(3.17) 



This condition corresponds to the Rayleigh criterium, i.e., different emission points cannot be resolved sufficiently to 
construct an "optical" image of source. Further, in most physical situations, L 3> a, so that 



\ki - /cal a < 1 



(3.18) 



This last inequality effectively resolves the artificial problem discussed at the end of the last section. The level spacing 
is sufficiently large that the analyzers do not detect particles coming from two different levels simultaneously. In other 
words, a measurement of the one-particle spectrum does not allow for a determination of the ca vity s ize. 

We now consider the two-particle spectrum starting from its Wigner operator representations ( ^.19[ ), or ( 2.20| ). An 
optimistic expectation is that the mean coordinates R = ±L will label the positions of the cavities, while coordinates 
z will be responsible for the internal dynamics of each of them, even after the operators Nk and Nki{Nk2 ~ (^feifea) 
have been averaged. We recall that the introduction of Wigner distributions depending on coordinates and momenta 
assumes that we have foregone a full quantum description of the particles. However, any quantum interference problem 
cannot even be formulated unless both final and initial states of the particles are described in terms of their quantum 
numbers. 

Both forms ( 2.19| ) and (2.2C) of the Wigner operator reproduce all four terms in ( [3.16 ), aft er av eragin g ove r the 



set o f qua ntum states of the two cavities. For example, consider Eq. (2.19). Substituting Eqs. (3.10) and (3.11 
Eq. ( 2.15 ), we obtain 



mto 



{NkANk2 ~ Sk,k2)) ^ 



E 

PI , ■ - ■ iP4 

Ni,...,N4 



/^t - \ ('^fei+'^Pl)(^fei+'^P2)(^fe2+'^P3)(^fe2+^P4) 

\"pi,Ari"p3,Af3"P4,JV4ap2,Af2/ ~ 



ATT^Aa'^ACLlki UJk2 ^/^pi^P2^P3^P4 



^ — ikiZi—ik2Z2 


^~ipi{Rl + ^} 




gV2(fl2 + ^) 




g-ip3(-R2-^) 




gip4(fli-^) 






Ni 




N2 




Na 





N4 
(3.19) 



The average of the product of Fock operators has four terms, given by Eq. ( 3.15| ), each corresponding to a specific 
combination of cavity states. The exponentials in ( 3.19| ) originate from the wave functions of the cavities, which set 
the limits of integration depending of the choice of cavity. For example, if Ni ^ N4 ^ —L and N2 =^ ^ +L then 
the limits of integration are 



L - a < Ri 



Zl 



< —L + a, L — a < Ri 



Zl 



< L + a: 



Z2 



-L - a < R2 — < -L 



Z2 



L-a<i?2 + — -<L 



(3.20) 



If Eq. (2.19) had represented only the first term of ( 2.34), then the first line of inequalities above would have contained 
only the cavity coordinate —L, while the second line would have contained only +L. This is not the case, and we 
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see that the coordinates Rj themselves fail to enumerate the cavities, contrary to our naive assumption. Thus, the 
requirement that the initial state be defined in terms of quantum states makes the Wigner distribution, which was 
introduced in a formal manner in Sections 2.1 and 2.2, a n inconvenient tool for handling in terfer ometry problems. 

If we consider all terms in the sum over Ni in Eq. ( 3.19 ), we recover the two-term formula ( 3.16 ) which was obtained 
from the original representation of the two-particle spectrum in terms of the locally defined quantum states. 



3.3 Hydrodynamics and interferometry 

The generalization of the two-cavity model to the freeze-out of an extended hydrodynamic system is evident: we 
consider a continuous set of decaying cells, taking into account the curvature of the freeze-out surface T{x) = Tc and 
the Doppler shift of the thermal spectrum of moving cells |^,^ . In addition, we neglected the effects of the multiple 
local emission, i.e., the possibility that two particles are emitted from the same elementary fluid cell. For the one- 
and two-particle spectra we then have [f7| 



,0^ 

dk 



J{k,k) , 



(3.21) 



and 



dNo 



J(fci,fci) J(fc2,fc2) + Re 



dkidk2 

respectively, where the emission function J(fci,A;2) is given by 



J{ki,k2) J{k2,ki) 



J{ki,k2 



dE^rr) n{k, ■ u{x)) e-('=-'=^)- 



(3.22) 



(3.23) 



In Ref. S, these equations were used to study interferometry for several types of one-dimensional flow. 



4. One— dimensional interferometry 



We shall analyze in this Section only the scale-invariant one-dimensional hydrodynamic regime, and neglect the 
boundary effects which were examined in Ref. In the central rapidity region, the scale-invariant solution does not 
differ much from the Landau model, provided that freeze-out takes place after sufficiently long evolution, so that the 
initial longitudinal size of the Lorentz-contracted nuclei in the c.m.s. is negligible. 

First, some notation: The parameters of the model are the critical temperature, Tc m-^), and the space-like 
freeze-out hypersurface, defined by — cc^ = = const. The rapidity y of a fiuid cell is restricted to ±Y in the 
c.m.s. We assume a Gaussian transverse distribution of hot matter in a pipe with the effective area S±_ nR\. The 
particles are described by their momenta 

fcf = {k°,kl,pi) = {mi cosh e. I, mjSinh6lj, pi) (4.1) 

where pi is the transverse momentum, Oi the particle rapidity in -direction, and mf — rn^ + pf \s the transverse 
mass. Let 



2a = 6li - 6*2 , 29 = + 02 (4.2) 

be the difference and the sum of particle rapidity, and q± = Pi — P2- The one-particle distribution has the form of a 
Bose thermal distribution of pions at temperature T ^ Tc. Since the m± values of interest are larger than Tc we may 
appro ximate the Bose distribution by a Boltzmann form, and use the saddle point method to estimate the integrals 



( 3.23 ). For the one-particle distribution these two steps yield 

ft! T.S I 777,1 / 

d0 



ft r5xmi r dy cosh(0i - j/) 6-"^ ^^'^(^i-^)/^^ ft rS^mi \[^^ e-"^'/^^ . (4.3) 
hdpi J_Y V rrii 



The general expression for J{ki, /c2) is 
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(toi +TO2) cosh(^^ — y) cosha — (mi — sinh(y — 0) sinha 



X exp 



1 

7; 



(mi + iF{nii~m2)) cosh{y — 9) cosha — (mi + i_F(mi + m2)) sinh(y- 



sinha 



(4.4) 
(4.5) 



where F = tTc- After some tedious algebra, using the saddle-point approximation, we obtain 



R{ki,k2) = Re 
1 



J{ki,k2)J{k2,ki) 
2TTfiT^S^Tc e 



g(^)g(|) 



X exp <; - — 



h{z) cos^^ + h{^)cos^^ 



X COS 



M.)sin^ + M>in^ 



H{z) + H{- 



+ G{z) + G{-) 



(4.6) 



where /i = (mim2)^/^ and z ~ (mi/m2)"'^/^, and we have introduced the functions 

n 2 



/i(z) 
tani?(z) 

5(^) 
tanG(z) 



4F^ sinh-" 



4F2(z2 -cosh2a)^ 



1/4 



2F(cosh2a - z^) 



^ Z 



1^2 



4F2 sinh^ a 



(z^ +cosh2a)2 + F2(z2 



1/2 



+ cosh 2a 



(4.7) 



The formulae (^^) and (4.7) are useful for a measurement of the longitudinal size in the case that the pions have 
unequal transverse momenta; the case of equal transverse masses mi = m2 = to_l was derived in Ref. The main 
result in that study was that the full longitudinal size of the freeze-out domain is not seen in the correlation function. 
Since 1-d expansion typically has large velocity gradients, the local thermal spectra of the slices with different rapidities 
do not overlap. The correlator measures the effective size of the fluid slice which forms the observable spectrum at a 
given rapidity 9 (or longitudinal momentum fc") 0: 



Ay = A — , 

m_L 



Axii 



cosh^ 



(4.8) 



These results seem to be in agreement with recent data [|l|,| , where the m^-dependence of the effective longitudinal 



size was examined. The physical origin of the dependence (4.8) is very simple: the larger the transverse momentum 



of a particle in the fluid, the more it is frozen into collective longitudinal motion, and the less the spectral pattern is 
broadened by thermal motion. Later we shall discuss what kind of parameters replace mj^ in the case of transverse 
hydrodynamic motion. 



The details of the correlator (3.22), in particular whether or not oscillatory behavior around unity may occur , hav e 
caused some controversy |lj,|l^. The co rrelat or derived in the formalism of the Wigner distributions (see Eq. ( ^.34[ )) 
does not reveal this behavior, while Eq. ( p. 22 ) permits such oscillations. We have already discussed why a naive usage 
of the Wigner formalism may distort interference effects. There is a simple argument why oscillations arc unavoidable 
in this class of hydrodynamic models: mathematically, only very restrictive conditions on the distribution functions 
will yield a positive definite result. Physically, when the velocity gradients are large (corresponding to small r in l~d 
expansion), two fluid cells with large relative velocity may still be close in coordinate space, even though the bulk of 
their spectra are Doppler-separated. We then effectively have two point-like sources, individually not resolvable, and 
the oscillations in the correlator are of the same origin as in HBT interferometry of a double stellar source. 

Up to now, we have assumed an instantaneous common proper time of emission for all particle types. However, it 
has been conjectured ||l3| that the hot expanding system may undergo a so-called "strangeness distillation," i.e. a 
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premature escape of the kaons from the hadronic gas because the K — tt cross-section is relatively small. We now 
discuss this phenomenon, and study how finite emission intervals may influence the interferometric data. 

1 /2 

A picture of sharp common freeze-out for all particles leads to the conclusion that 7 = a^ffirn^ is independent of 
TO^, and that 7^ = ^k- This is indeed observed experimentally 1^, suggesting that the kaon fluid remains coupled to 
the pions right until their common (sharp) freeze-out. If we allow for a dynamical decoupling of the kaons from the 
pionic fluid, then the full kaon spectra will result from a gradual emission of the kaons over some interval ro < r < r^r. 
For simplicity, we shall assume that the emission function is given by the unweighted average 

{J)r = rdrJir) , (4.9) 

- To Jrn 



where J is calculated numerically from Eq. ( 3.23 ), and is the pion freeze-out time. In the integral (4.£), we 
assume that the temperature is a function of r, according to tT^ = const. In subsequent calculations we shall take 
= 130 MeV at = 30 fm. 

In Fig. |l|, we show the correlator as a function of Ak^^ /m^/'^. The various curves correspond to different values of 
m± . The curves marked "s" are the result if we take tq = for the emission function. They coincide almost perfectly 
with each other, in agreement with Eq. ( |4.8| ). The curves marked "g" correspond to a gradual emission with tq ~ 7 fm 
(To = 180 MeV). Clearly, the m^^^-scaling is violated for this type of emission. NA44 data gives no indication of 
such scaling violation. A systematic and parallel study of pion and kaon interferometric source sizes as a function of 
m± will provide a sensitive test for strangeness distillation in hadronic matter at RHIC. 



5. Two— and three— dimensional interferometry 

In this Section, we discuss the interferometric measurements for two- and three-dimensional hydrodynamic flow. 
To motivate this discussion, we begin by recalling that the most important lesson from the case of a one-dimensionally 
expanding source is that there is no universal prescription for decoding interferometric data. Thus, in interpreting 
this data, we must first determine which model is applicable to the data. The model should depend on a minimal 
number of parameters, and exhibit some specific and clear behavior that is qualitatively recognizable in the data. 
Only then may we hope to extract these unknown parameters from interferometry. For example, in 1-d hydrodynamic 
motion of the Bjorken or Landau type, the signatures are the plateau in the central rapidity region, an increase of the 
correlator width with increasing particle pair momentum, independence of the correlator width upon the sum of the 
rapidities, independence of the product aeffy/mj_ on the mj_, and an increase of the correlation function at rapidities 
that exceeding the plateau width. If such behavior is found, one may use the corresponding formulae and appropriate 
variables in fitting the data, and, using the full dynamical description of the evolution, determine parameters of the 
expanding system. An example of such a strategy was described in Ref. The abundance and clarity of signatures 
for a boost-invariant-like scenario are the result of the homogeneity of the flow parameters along the expanding hot 
pipe. The contributions from all the slices add to intensify the signal. Mathematically, this type of flow allows for 
rather simple analytic formulae, and since we are free to choose sufficiently large m±, the saddle point approximation 
is very accurate. This fact is important, because only with analytic formulae can one obtain an understanding of the 
dependence of the correlator on the physical parameters. 

Unfortunately, this attractive feature of the boost-invariant solution for one-dimensional flow is absent in two- 
and three-dimensional flows, even if the radial motion is very strong. For radial flow, the angular dependence of 
the Cartesian velocity components is weak, and the localization of the momentum spectrum is far less than in one- 
dimensional expansion. These small gradients mean that a saddle point integration will be a bad approximation, and, 
consequently, we have no simple dependence of the correlator upon the flow parameters. While even moderate radial 
flow does obscure the true transverse source size, the spectrum is not sufficiently localized to allow one to obtain a 
simple formula for the correlator (at least in terms of standard variables) . 

Let us illustrate this using two unphysical, but exactly calculable examples. The first is a scale-invariant two- 
dimensional expansion in the absence of longitudinal flow, as in the explosion of a long, thin fllament. The second 
example is a scale-invariant three-dimensional expansion, corresponding to a point-like explosion. An advantage of 
these models is that the hydrodynamic equations are exactly solvable, and we are guaranteed dynamical consistency 
between the velocity field and the shape of freeze-out surface. This is an important requirement for any model 
calculation. 



5.1 The explosion of a long filament 
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We consider the transverse expansion of a filament of length L, much larger than the transverse freeze-out radius. 
For this case of purely cylindrical expansion, a convenient parametrization of the coordinates is 



= (t cosh/3, T sinh/3 cos?/", r sinh/3 simp, z) 
while the temperature and the velocity fields may be written as 

= const 

— (cosh/3, sinh/3 costA, sinh/3 s\\\ip, 0) , 



(5.1) 



(5.2) 
(5.3) 



where (3 turns out to be the radial rapidity of the fluid element. An element of the freeze-out hypersurface is given 



by dT.^' 



i'^ sinh (3 dz d(3 dip. 



The emission function J{ki, fe) may be expressed using the integral 



J 



2tt 



sinh P d/3 dip exp 



a cosh /3 + &i sinh/3 cos-0 + &2 sinh /3 sin 



27r 



sinh/3d/3e~''™"^'^/o(Vfe? + 6i sinh/3) 



27r 



hi 



and derivatives of with respect to its parameters. Introducing the shorthand notation 



H = 



1 + 



we find 



J{ki,k2) 



T 



{m,H/Tf 



T 



2to? 



J(fc2,fcl) 



which immediately gives us the one-particle spectrum of the model: 

dk [m,JT) 



T 



1 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



Similar to one-dimensional boost invariant flow, we obtain a plateau for cylindrical boost invariant expansion, but 
now in the radial rapidity distribution. The dependence of the spectral density on to^ is noteworthy: the localization 
of the spectrum due to radial flow is more pronounced for greater m^. Particles with m^/T 3> 1 are strongly frozen 
into collective radial flow. For this reason we may expect the width of the correlator in the transverse direction to be 
defined by m^, rather than m^. 

The expression for the normalized two-particle spectrum is 



C(fci,/c2) - 1 = Rc 



'2m^{H~l)/T 
IP' 



m,H/T+l 



n 2 r 



mjT+1 



1 



8m? 



(5.8) 



Its explicit real form is of little practical value, except to estimate the reliability of the saddle-point approxim ation 
for the two- and three-dimensional flows. We shall return to this point later. Instead, let us first approximate (5.S) 
in the case of small differences in radial rapidities of the particles. For the sideways and outward directions we obtain: 



C(Q) - 1 



cos 



Tt 
2m1 



Qout 



for Qside = 0, and 



C(Q) - 1 = e-^^'^"-/'^ 



cos 



Tt 



—CP' 

^ side 

for Qout = 0. We then obtain the same radii for the sideways and outward directions: 



■exp 



(5.9) 



(5.10) 



(5.11) 
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From symmetry, it is not unexpected that the radii should be the same for sideways and outward directions. The radius 
Rexp is defined by the exponential, and dominates the shape of the correlat or fo r Tt > 1, while Rcos, defined by the 
zero of the cosine function, dominates for Tt < 1. Formally, the expressions ( ^.9[ ) and ( 5.10| ) are valid at small relative 



transverse rapidities of the pions. They very much resemble the basic formula for longitudinal expansion. However, 
one should keep in mind that the latter are obtained wit hou t the saddle-point approximation. This approximation 



fails in the cylindrical case because the integrand of Eq. (5^) varies only slowly with the azimuthal angle -0 over the 



entire domain of the integration. To clarify this statement, we have estimated numerically by how much the integral 
is affected if we retain only a small ■(/'-interval about the maximum of the integrand. Even for unreasonably large 
ratios m^/T ^ 10, we must integrate over a large interval 1-01 < tt/S in order to obtain 90% of the exact value of 
dN/dk. 

5.2 The explosion of a point-like source 

We next consider the case of spherical expansion, and parameterize the coordinates as 

x^ = T (cosh/3, sinh/5 sin 6* cosip, sinh/3 sin 6* sin0, sinh/3 cos6) , (5.12) 
and the temperature, velocity field, and freeze-out hypersurface volume element can be written as 

tT = const, u'' = x^/r , dW = u^" sai\? (3 saiO d(3 dO d^p . (5.13) 
Once again, to calculate the emission function J(A:i, ^2), we use the auxiliary integral 

/>oo />7r />27r 

J= Sinh^ (3d(3 sin 9d0 ^^g-acosh/3+h3sinh/3cose+si„h/3sine(6i cosV>+b2 sinV') ^ 



Jo Jo 



An 



K,U a^-til-bl-hl) , (5.14) 



were Ki is the modified Bessel function. It is now straightforward to find the one- and two-particle spectra: 



k^'^ ^ Ji^k^k) ^ Anr^TK^C^) , (5.15) 
dk ^ 



C(fci,/s2) - 1 = Re 



K2{H) 



H^K2{m/T) 



1 ^ 



t2 



We see that spherically symmetric boost-invariant expansion washes out any inhomogeneity of the local thermal 
spectrum. This is a distinctive feature of the spherical model. Approximating the correlator in the case of small 
differences in radial rapidity of the particles, we obtain 

0,0,-1 = e--«=.-,i±g^eo.rg,5Il + I^A . ,=.16) 



l-H5r2T2Q2/2m2 Vto2' 2 2 
The radius which may be extracted from this correlator is 



Rexp = T\ — , (5.17) 
" TO 



and the model has no large parameter which could allow an asymptotic estimate of the emission function. 
6. Conclusion 

Preliminary results recently reported by the NA35 and NA44 collaborations seem to confirm the main predictions of 
interferometry for a hydrodynamic one-dimensional flow scenario of the heavy ion collision. The intensity correlations 
provide a multidimensional test of the source, and the fact that so many parameters coincide can hardly be accidental. 
Thus, we have strong evidence that even at SPS energies (s^/2 ~ 10 A • GeV), a hydrodynamic regime develops, and 
that the freeze-out takes place during a short interval. This collective behavior can be expected to occur at RHIC 
energies, and thus it is highly desirable to continue developing the formalism. 
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In this paper, we have emphasized the most important physical aspects of the theory, analyzed the reliability of its 
different modifications known in the literature, and tried to isolate some controversial points. We argue that, strictly 
speaking, interferometry does not permit the initial data to be given in terms of semi-classical distributions, unless 
this description is augmented by a clear indication of the length scale that defines the quantum states of the particles. 
Our conclusion here is that the traditional operator approach, based on the precise definition of the particle states, 
prov ides a firm footing for the calc ulation of the two-particle spectra. In this case, distributions such as n{p, N) in 
Eq. ( |3.15D , and n{k, x) in Eq. ( |3.23| ) can be thought of as on-shell Wigner phase-space distributions of the preceding 



kinetic stage. For any type of hydrodynamic expansion, no principal or technical problems are anticipated. 

We have extended previous calculations ||] for the longitudinally expanding system to the case of unequal transverse 
momenta, and have shown that the HBT correlator may carry information about strangeness distillation. From the 
NA44 data analysis, it seems that such distillation does not occur at SPS energies. We further considered transversally 
expanding sources, and demonstrated that in this case one has to change the parameterization of the correlator. 

We emphasize that, while HBT for hydrodynamic sources is well understood as a physical phenomenon, it is always 
a problem to choose the "correct" model for fitting the data. We impose the condition that the model should allow 
one to recognize it via a qualitative analysis. Only then can one hope to understand what parameters are responsible 
for the correlator behavior, and find their values by fitting the data. Only after such an analysis has been performed, 
does it become possible to use the dynamical equations of the model, and trace the freeze-out parameters back to 
the earlier stages of evolution and, eventually, to an estimate of the energy density. Practically, this requirement 
means that we must begin with an analytic expression for a solution of the relativistic hydrodynamic equations. This 
guarantees the consistency between the shape of the freeze-out surface and the velocity field. Unfortunately, analytic 
solutions for the case of three-dimensional expansion are not yet known (the equations are nonlinear, and the initial 
data is not well known). A reasonable analytic approximation will do, but to our knowledge, a suitable expression 
has not yet been derived. An approximate formula describing a realistic, expanding system at the freeze-out stage is 
an important problem for boson interferometry at RHIC. 
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Appendix 

Descriptions of nuclear collisions in terms of semi-classical kinetic equations rely heavily on the classical nature of 
the source that emits the final-state particles. Such sources may be expressed naturally in terms of Wigner functions, 
and it is therefore of interest to what extent these functions may be used directly to solve the interferometry problem. 



To our knowledge, equations like (2.34) and (2.35) were originally derived in a model with classical sources. Here we 
revisit this derivation. 

Let ifiix) be a quantum pion field that is emitted by a classical current j{x). The corresponding solution of the 
Schrodinger equation for the state vector is a coherent state |<I>), with the property that 

m0{xm - v,i{x) , (A.i) 

where the classical field ipci{x) obeys the inhomogeneous wave equation with source j{x): 

{a + m')^,i{x) = j{x) . (A.2) 

For space-time regions outside of where the currents j{x) are localized, (pdix) obeys the homogeneous wave equation, 
and may be expanded in the plane-wave modes of Eq. (|]^) : 

Vciix) - J d^ka^J^ix) . (A.3) 

The coherent state |<I>) may be obtained by solving the following equation for each mode: 

and, after normalizing to = 1, one obtains the representation 
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E 



\nk) 



i"£r/2 



ni,=0 



From the equation of motion for the classical field, ( [A.2| ), one immediately obtains 
so that the coherent state |<i>) may be rewritten as 



1$) 



exp 



^3, ljK:fc)l' 



exp 



Ai 



(27r)3/2(2cjfe)i/2 



|0) 



(A.5) 



(A.6) 



(A.7) 



We may also write |<i>) = >5'|0), and it can be shown explicitly that the expression ( |A.7| ) emerges if the S'-matrix is 
taken in the form of a normal product. 

Evaluating the one- and two-particle spectra is now a question of determining the right hand sides of 



^(^1-^2)))) 



($|it it ir i.- I*) 

^ ' fci fc2 ^=2 fcll ' 



(A.8) 
(A.9) 



where ((...)) denotes an ave rage oyer the statistical ensemble of classical currents. These equations may be obtained 
immediately from Eqs.(2.4) and (2.5), if the density matrix is taken of the form 



(A.IO) 



Since the coherent state |<I>) is the eigenstate of every annihilation operator A^, with eigenvalue ag, we immediately 
obtain the final result 



((iV,-)) = {{ataj:)) = (( 
{{Nj:^iN^,^-S{h-k,)))) 



Qfi ttr oir Cir 

fci k-2 ki k2 



j*{ ujk,k)j{ujk,k) 

II 3 * (^/ci , ki )j *{ujk2, k2 )i {^kt ,ki)j {ujk2 , k2 ) 



The Fourier transform of the currents is defined in a standard way: 

j{k) = j{ko,k) = J d^xj{x)e^^^ , 
and we therefore obtain our previous answer for the one-particle spectrum 



The two-particle spectrum now reads 



2a;fe(2^)3 



-i k-{x~-y) 



{{Nj:^{Nf:^ - 6{k, - fc^)))) = / '^^'^^^''1'/ e-'=-(-^)-'=-(-'-^') {{f{x)f{x')j{y)j{y')}} 



(A.ll) 
(A.12) 

(A.13) 

(A.14) 

(A.15) 



i.e., we recover the well known result that a coherent source does not lead to any quantum interference effects. In 
other words, quantum evolution follows a single trajectory, which ends up as a pure coherent state of the pion field. 
After we have averaged over this state, further consideration of the two-pion amplitudes and their interference is 
impossible. However, we note that there remains the opportunity to account for statistical effects inherent in the 
distribution of the classical currents. An example is given in Ref. coherent light propagating in a fluid induces 
(classical) polarization currents in the molecules. The latter become sources for a secondary (scattered) field. Intensity 
correlations of this field allow one to study the multi-particle distributions in the fluid. For simple fluids, the exact 
solution of the inverse problem is possible 0| . 

One may formally consider the S'-matrix in Eq.( [A.7| ) as the perturbation series (assuming, e.g., that the currents 
are weak). For such an expansion it is convenient to use the S'-matrix in the form of a T-ordered exponent: 
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S = Texpl -i / d^x[f{x)ip{x) + j{x)ip'' {x)]} . (A.16) 



If only terms up to second order in the currents are retained, then the coherent state is effectively replaced by a 
two-particle state. This state is no longer the eigenstate of the annihilation operator, and the matrix element that 
must be evaluated now reads 

d^xd'yd'x'dy {{3*{x)r{x')j[y)j{v'))) (0|^(x)^(a:')it/t/4i£^(pny)^Hy')|0) (A.17) 

Commuting the creation and annihilation operators with the field operators we arrive at 



/ 



\ ikix+ik2x' , ikix' +ik2X^ ikiy+ik^v' _i_ pikiy' +ik2y] 

d\d'yd'x'dy{{3*[x)3*{x')j{y)3{v'))) ^ ^^1^ — ^ . (A.18) 



Choosing instead of {x^x') the coordinates (i?, z), we obtain the form 

4fc0fc°^^ = / d^R^d^R 
dki dky J 



dkidk2 

where 



F{R,,k,)F{R2,k2) + e-(«-^^)-('^--^^)F(i?i,^^l±^)F(i?2,^^) 



(A.19) 



F{R,P) = J -^e-^P- j*{R + z/2)j{R~z/2) (A.20) 

is the Wigner representation of the current product. In calculating the average ((...)), one usually assumes that the 
product of four currents factorizes into the binary products. 

Clearly, this way of reasoning contradicts the concept of the classical current itself: if the field is classical, then 
the number of quanta is undefined, in principle. The selection of fluctuations with the exactly two emitted quanta is 
an additional measurement of the intermediate state of the system. This measurement prepares a new state of the 
system and removes all information about the classical nature of the emitting currents. The further evolution can 



now follow two different trajectories. This is exactly why Eq. (A. IS) has acquired a typical interference structure, 
even before the particles histories were traced back to their sources. 

On the other hand, if we can select a state with two quanta in the emission field, then the current is essentially 
a quantum mechanical transition current. In this case we must take into account an explicit backward reaction in 
the emitting system, and we arrive naturally at the picture which has been already considered in connection with 
photon emission. We conclude that there is no consistent way to express the correlations in terms of off-shell Wigner 
functions. 

* E-mail: makhlin, gene, and welke@rhic.physics.wayne.edu, respectively. 



[1] The NA35 collaboration: R. Morse et ai, Preprint LBL-36062 (1994). 

[2] The NA44 collaboration: H. Baker et ai, Preprint CERN-PPE/94-119 (1994). 

[3] For an overview, see, for example: M.I. Podgoretskii Sov. J. Part. Nucl. 20 (1989) 266. 

[4] A. Makhlin, Sov. Phys. JETP 35 (1972) 478. 

[5] S. Pratt, Phys. Rev. Lett. 53 (1984) 1219. 

[6] G. M. Welke, H. W. Barz, G. F. Bertsch, P. Danielewicz, and H. Schulz, in Proc. of the 10th Winter Workshop on Nuclear 
Dynamics, January 15-22, 1994, Snowbird, Utah (World Scientific, Eds. J. Harris, A. Mingerey, and W. Bauer) p. 93. 

[7] A. Makhlin, and Yu. Sinykov, Sov. J. Nucl. Phys. 46 (1987) 354. 

[8] V. Averchenkov, A. Makhlin, and Yu. Sinykov, Sov. J. Nucl. Phys. 46 (1987) 905. 

[9] K. Kolehmainen, and M. Gyulassy, Phys. Lett. B 180 (1986) 203. 
[10] M. Gyulassy, S.K. Kaufrniann,and L.W. Wilson, Phys. Rev. C 20 (1979) 2267. 
[11] A. Makhlin, JETP Lett. 46 (1987) 55. 
[12] A. Makhlin, Sov. J. Nucl. Phys. 49 (1989) 238. 
[13] D.K. Srivastava, and J.I. Kapusta, Phys. Rev. C 48 (1993) 1335. 



17 



[14] H. Sorgo, H. Stockor, W. Grcincr, Nucl. Phys., A498 (1989) 567c. 
[15] J. Rafelski, Phys. Reps. 331 (1982). 

[16] See references in W.A. Zajc, in Particle Production in Highly Excited Matter," H. Gutbrod, and J. Rafelski, eds., NATO 

ASI series B303 (Plenum Press, 1993), p. 435. 
[17] M. Born and E. Wolf, Principles of Optics (Pergamon, Oxford, 1964). 

[18] A. Timmermann, M. Plumer, L. Rasumov and R.M. Weiner, Preprint DMR-THEP-93-8/W, <hep-ph-9405232> (1994). 
[19] S. Chapman and U. Heinz, <hep-ph-9407405> (1994). 



Figure Caption 

FIG. 1. The correlation function C2, as a function of the rcscalod longitudinal momentum difference, for different values of 
m±. The curves marked "s" correspond to sharp emission, and those marked "g" are for gradual emission. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9412359vl 
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